Effects of polymer polydispersity on the phase behaviour of colloid-polymer mixtures 
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We study the equilibrium behaviour of a mixture of monodisperse hard sphere colloids and poly- 
disperse non-adsorbing polymers at their #-point, using the Asakura-Oosawa model treated within 
the free-volume approximation. Our focus is the experimentally relevant scenario where the distri- 
bution of polymer chain lengths across the system is fixed. Phase diagrams are calculated using the 
moment free energy method, and we show that the mean polymer size £ c at which gas-liquid phase 
separation first occurs decreases with increasing polymer polydispersity 5. Correspondingly, at fixed 
mean polymer size, polydispersity favours gas-liquid coexistence but delays the onset of fluid-solid 
separation. On the other hand, we find that systems with different S but the same mass-averaged 
polymer chain length have nearly polydispersity-independent phase diagrams. We conclude with a 
comparison to previous calculations for a semi-grandcanonical scenario, where the polymer chemical 
potentials are imposed; there it was found that fluid-solid coexistence was favoured over gas-liquid 
in some areas of the phase diagram. Our results show that this somewhat counter-intuitive result 
arose because the actual polymer size distribution in the system is shifted to smaller sizes relative 
to the polymer reservoir distribution. 

PACS numbers: 



I. INTRODUCTION 



Colloid-polymer mixtures have a wide variety of practi- 
cal uses: polymers can be added in low concentrations to 
colloidal suspensions to modify their properties 0, 0, |j| , 
or conversely colloids can be added to polymer melts as 
"fillers" 4] or to polymer solutions to produce gels [f| . In 
the former case, it is well known that the presence of the 
polymer can induce an effective colloid-colloid depletion 
attraction 0, B 0- Unlike in atomic systems, both 
the strength and range of this interaction can be tuned 
by varying polymer concentration and size. The implica- 
tions of this for the phase behaviour of the mixture are 
of fundamental theoretical interest, and have been the 
subject of much research. 

We focus on the simplest situation in this paper, where 
the colloids have hard interactions with each other and 
with the polymers; this requires in particular that the 
polymers be non-adsorbing on the surface of the colloids. 
We further restrict ourselves to 0-point conditions, where 
polymer-polymer interactions can be neglected to a first 
approximation. When the polymers are modelled as ran- 
dom walks on a lattice, this scenario is amenable to direct 
numerical simulation |lCI However, most work has 
focussed on a further simplification, proposed by Asakura 
and Oosawa (AO) [a, ancl Vrij Q , where the polymers 
are replaced by spheres with radius equal to the radius 
of gyration of the original polymer chains. These "poly- 
mer spheres" are then assumed to interpenetrate freely 
with each other while having a hard excluded volume 
interaction with the colloids. This AO model has been 
studied in great detail, by thermodynamic perturbation 
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theory 0,EJ0> free- volume theory 0] (see below), 
density functional theory 0, 0, , and Monte Carlo 
simulations 0, 0, |2jj. Note that the naive choice of 
setting the radius of the effective polymer sphere equal 
to the polymer radius of gyration can be improved upon; 
we will return to this point in the discussion in Sec. IIVI 
In addition to the simplifications already discussed, the 
AO model ignores the fact that there is, in real colloid- 
polymer mixtures, invariably an (essentially continuous) 
spread of polymer chain lengths. The effects of such poly- 
mer size polydispersity on the phase behaviour have been 
studied by some researchers. As described in more de- 
tail below, however, the experimentally relevant situa- 
tion where the full polymer size distribution in the sys- 
tem is conserved during phase separation remains to be 
understood. This is the issue which we address in this 
paper. Specifically, we study the AO model within the 
free- volume approximation of [l5j , for the case where the 
colloids are monodisperse, with diameter <r c , while the 
polymer spheres are polydisperse with density distribu- 
tion Pp((Tp) = p p n p (a p ). Here p p is the polymer number 
density and n p (o~ p ) is the normalised polymer diameter 
distribution. We take the colloid diameters a c as our 
unit of length so that the polymer-to-colloid diameter ra- 
tio £ coincides with the dimensionless polymer diameter, 
£ = fjp. All densities are made dimensionless by multi- 
plying by the volume of a colloid particle, and quantities 
with the dimension of energy are measured in units of 
k B T. 

For the case of monodisperse polymers, it is well 
known that the phase behaviour depends stron gly on the 
polymer-to-colloid diameter ratio (see e.g. [TEl l2l|). For 
small polymers only coexistence between colloidal fluid 
and solid phases is observed. (The term "fluid" is used 
here because there are no distinguishable gas and liquid 
phases in such a system.) For larger polymers, on the 
other hand, with £ above some threshold value £ c , a re- 
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gion of gas-liquid phase separation can occur in the phase 
diagram. When generalising these considerations to the 
case of polydisperse polymers, there is no longer a single 
polymer size but rather a continuum of different diame- 
ters <7 p . A choice therefore needs to be made on how to 
define appropriately the typical polymer size that deter- 
mines the phase behaviour. Ideally, one would like this 
definition to be such that phase diagrams become nearly 
independent of polymer polydispersity and can therefore 
be inferred from the corresponding monodisperse refer- 
ence system [22^. 

A naive definition of the typical polymer size is the 
mean polymer diameter, £g = J da p a p n p (a p ) . This is, 
in fact, independent of the diameter polydispersity S, de- 
fined as the standard deviation of n p (cr p ) normalised by 
its mean. We will see below, however, that this provides 
a poor mapping from polydisperse to monodisperse sys- 
tems, with phase behaviour at fixed £g strongly depen- 
dent on S. Following [22j, we have therefore investigated 
other definitions of typical polymer size. These include 
the number average £jv 5 defined by 



J da p a1n p {a p ) 



(1) 



This is proportional to the average polymer chain length, 
which scales as a p . The mass average £m is also an av- 
erage chain length, but with the contribution from each 
polymer size weighted by an additional factor of chain 
length (i.e. mass): 



f 2 



J da p a^n p {a p ) 
f da p a%n p (o- p ) 



(2) 



A final definition, the virial average £b, is derived from 
the second virial coefficient for the polymer-colloid inter- 
action. A polymer's centre of mass is excluded from a 
sphere of diameter 1 + tr p around a colloid of size c c = 1, 
and £b is defined as the monodisperse polymer diameter 
giving the same average excluded volume, viz. 



(1 



S,b) 3 = J da p (l + a p ) 3 n p {a p ) 



(3) 



In a monodisperse system all of these definitions are of 
course equivalent, so that £s = £jv = £m = £b- In 
a polydisperse system they respond differently to the 
shape of the polymer diameter distribution, and we will 
try to identify which definition is best suited to charac- 
terising the phase behaviour in a manner that depends 
only weakly on the polymer polydispersity. In assessing 
this dependence we also have a choice for an appropri- 
ate measure of polymer concentration. This could be 
e.g. the total polymer number density p p = J da p p P (c p ) 
or the effective volume fraction of the polymer spheres, 
4> p = Jda p a p p p (a p ); values of <fi p of order unity cor- 
respond to the crossover from the dilute to the semi- 
dilute polymer regime. For a fixed diameter distribution 
n p ((Tp), p p and <p p are proportional to each other, but 
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FIG. 1: Phase diagram sketch of the phase behaviour for a 
semi-grandcanonical colloid-polymer mixture, following [2^1: 
F: fluid, G: gas, L: liquid, S: solid. The two graphs show 
the colloid volume fraction on the horizontal axis and on the 
vertical axis the polymer reservoir density pp (left) and the 
polymer density p p in the system itself (right). The dashed 
and solid lines correspond to weak (z — > oo, 5 = 0) and strong 
(z — 5) polymer polydispersity, respectively; the virial aver- 
age polymer size is £b = 0.4 in both cases. The triple points 
are marked by horizontal dotted lines on the left; they are 
mapped to three-phase GLS triangles on the right. 



when comparing distributions with different shapes the 
difference in the proportionality coefficients matters. We 
choose mostly to work with </> p , following Warren's sug- 
gestion that this produces phase diagrams with a 
relatively weak dependence on the polymer diameter dis- 
tribution. This will allow us to assess how far Warren's 
insight generalises beyond the binary mixtures of small 
polymers which he considered. 

As far as previous work on polymer polydispersity ef- 
fects in colloid-polymer mixtures goes, we have already 
mentioned Warren's work for mixtures of polymers of 
two different sizes j2^|. This explicitly accounted for the 
fact that the overall polymer density distribution in the 
system - which in the polydisperse case we will write 
as p p (<T P ) - is imposed by the experimental conditions. 
Studies of polymers with continuous size polydispersity, 
on the other hand, have been restricted to the semi- 
grandcanonical case |23ll24j. (An exception is very recent 
work |2^| on interacting polymers; see Sec. ED Here the 
system is thought of as connected to a large polymer 
reservoir which fixes the chemical potential p p (o~ p ) for 
each polymer size a p . The polymer density in the sys- 
tem is determined only indirectly, and has to be deduced 
a posteriori. Within the AO approach the polymers are 
treated as ideal so that the reservoir density distribution 
and chemical potentials are related by p p (o~ p ) = e Mp (°p'. 
The decomposition of the chemical potential into ideal 
and excess parts, p, p {a p ) = In p p ((7 p )-|-p; p x (<7 p ) then shows 
that the polymer density distribution in any of the phases 
within the system itself can be written as 



Pp 



(4) 



Here the /Ltp X ((T p ) are the excess polymer chemical poten- 
tials of the phase being considered. 

Fig-dshows a sketch of the results obtained in Ref. [24| 
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for such a semi-grandcanonical scenario, with polymer 
reservoir density distributions of Schulz form (see Sec. 11111 
below). The graphs compare the behaviour for weak 
(z — > oo, i.e. 5 — 0, dashed lines) and strong (z = 5, 
solid lines) polydispersity. The means of the Schulz dis- 
tributions were chosen so that the virial average poly- 
mer size £b = 0.4 of the normalised reservoir size dis- 
tribution was the same in both cases. On the left, the 
polymer reservoir density p p is sketched against the col- 
loid volume fraction. Using Q, these results can then 
be mapped onto a representation in terms of the actual 
polymer density p p in the system, as shown on the right. 
This corresponds to fixing the total number of polymers 
in the system while the polymer size distribution contin- 
ues to vary across the phase diagram, being controlled by 
the chemical potentials - or rather the chemical potential 
differences - imposed via the reservoir. 

In both representations Fig. H] shows that, as a general 
rule, increasing the polymer polydispersity 5 increases 
the size of the gas-liquid region in the phase diagram. 
As seen most clearly in the plot on the right, however, 
there are also regions in the phase diagram (at low p p ) 
where gas-liquid coexistence is suppressed with increas- 
ing 6, and fluid-solid phase separation is favoured. This 
seems counter-intuitive: size polydispersity in the col- 
loids, for example, is known to have the opposite effect, 
dela ying fluid-solid relative to gas-liquid phase separa- 
tion [2l|, |23, |23 ■ One of our aims will be to clarify this 
issue, and in particular the effect of the varying polymer 
size distribution in the semi-grandcanonical setting. 

The effects of polymer size polydispersity have also 
been analysed in terms of the (pairwise part of the) deple- 
tion interaction between colloids 28,29,30,31]. Concep- 
tually, the fixed polymer size distribution in these stud- 
ies means again that constant imposed polymer chemical 
potential are being considered; see e.g. Ref. 11511331 . Fi- 
nally, we mention the related work of Ref. [33. This 
also used a semi-grandcanonical framework, but start- 
ing from a different physical motivation: polymers were 
regarded as monodisperse with regard to chain length, 
and size polydispersity in the equivalent polymer spheres 
was considered as arising from compression of the poly- 
mers. In this case the semi-grandcanonical approach is 
the physically relevant one since polymer spheres of dif- 
ferent size can be transformed into each other. (A similar 
scenario arises when the polymers are replaced by worm- 
like micelles The effects of polymer compressibility 
on fluid-solid phase boundaries were found to be rather 
small up to size ratios £ « 0.5 [13 • We therefore ignore 
them here and focus on size polydispersity arising from 
a distribution of polymer chain lengths. The experimen- 
tal situation is then normally the canonical one, where 
the overall distribution of polymer sizes in the system is 
fixed. 

Below we describe briefly the free energy which we use 
to model our colloid-polymer mixture, and the numerical 
method by which phase diagrams are obtained fSec. ITTfl - 
Our results for the effects of polydispersity on the phase 



behaviour are given in Sec. IIIII separately for the ex- 
perimentally relevant case of fixed polymer size distribu- 
tion and for the semi-grandcanonical scenario where we 
compare with previous studies p4j ]. Conclusions and an 
outlook towards future work can be found in the final 
section. 



II. FREE ENERGY AND NUMERICAL 
METHOD 

In a recent paper [2l| we derived a free energy expres- 
sion for an AO system where both polymers and colloids 
are polydisperse, by extendin g th e free-volume approxi- 
mation of Lekkerkerker et al |l5j . This free energy can 
be written in the following form: 

/ = Jf + f? + fS + / da p Pp( CT pK>p) (5) 

The first two contributions are the ideal free energies of 
a mixture of polydisperse colloids and polymers, respec- 
tively. The next term, f%£, is the excess free energy of 
a system of pure (polydisperse) hard-sphere colloids. In 
the last term, which represents the colloid-polymer in- 
teraction, Pfe(o-p) is the excess chemical potential of a 
hard-sphere particle of diameter a p in the pure colloid 
system. 

Evaluation of the free energy (JSJ requires as input only 
the properties of the pure colloid system, i.e. its excess 
free energy f£*. We therefore need to assign appropriate 
expressions for in the colloidal fluid (or gas/liquid) 
and solid phases. For the fluid part of the excess free 
energy the most accurate approximation available is the 
BMCSL equation of state [35|, |3g while for the solid we 
adopt Bartlett's fit to simulation data for bidisperse hard 
sphere mixtures [13, H^- Our previous work [26l |2?| 
on polydisperse hard spheres has shown that with these 
free energy expressions quantitatively accurate fits to 
simulation data arc obtained. Both free energy expres- 
sions depend only on the moments p C i — J da c p c {a c )a l c 
(i = ... 3) of the colloid density distribution. As a 
consequence the excess chemical potentials become third 
order polynomials, 

X fex 

where the i are the excess moment chemical potentials 
of the pure hard-sphere system, = df^/dp C i. (As 
explained in Ref. [21| , we always evaluate these from the 
BMCSL free energy.) The interaction term in JSJ then 
simplifies to 

J do- p Pp(o-p)/4*((Tp) = Y^Vte,iPpi ( 7 ) 

i 

where the p P i = J da p p p (o~ p )a p (i = . . . 3) are the mo- 
ments of the polymer density distribution. 
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We will deal with the particular case where the colloids 
are assumed to be monodisperse. With the colloid diam- 
eter set to unity as assumed, all colloid density moments 
then become identical, p C i — p c . Writing out the ideal 
contributions explicitly, the free energy (JSJ thus takes the 
form 

/ = J da p p p (a p )[\n p p (a p ) - 1] +p c (lnp c - 1) 

+ /hs(Pc)+£^*>cK* (8) 



Here we have highlighted that (and therefore the 
Mhs i) now oru y depend on the colloid density p c , or equiv- 
alently the colloidal volume fraction <p c ; for monodisperse 
colloids and in our units the two quantities are identi- 
cal. The excess part of the free energy (JHJ thus only de- 
pends on the colloid density p c and the moments p P i of 
the polymer density distribution. As in [2l| for the con- 
verse situation of polydisperse colloids and monodisperse 
polymers, we can regard these quantities as moments of 
an enlarged density distribution (p c , p p (a p )). Because 
only a finite number (five) of such moments are involved, 
we thus have a truncatable free energy [39|. The phase 
equilibrium conditions associated with this free energy 
can then be solved using the moment free energy (MFE) 
method |H |M ES El El • The moment free energy al- 
lows one to map the full free energy JSJ, with its depen- 
dence on all details of p P (cr p ) through the ideal part, onto 
a MFE depending only on the moments p P i and p c . From 
the latter, phase behaviour can then be found by the 
conventional methods for finite mixtures, treating each 
of the p„ i as a number density of an appropriate "quasi- 
species" [39j. For truncatable free energies this locates 
exactly the cloud points, i.e. the onset of phase separation 
coming from a single phase, as well as the properties of 
the coexisting "shadow" phases that appear. Inside the 
coexistence region, one in principle needs to solve a set of 
highly coupled nonlinear equations, and the predictions 
derived from the MFE are only approximate. However, 
by retaining extra moments, increasingly accurate solu- 
tions can be obtained by iteration [3^, \43L EH |45| . Using 
these as initial points, we are then able to find the, for 
our free energy, exact solutions of the phase equilibrium 
equations. 

The MFE method is designed for the physically real- 
istic scenario where the overall polymer density distribu- 
tion is conserved when the system separates into two or 
more phases. However, it can also be exploited to analyse 
semi-grandcanonical scenarios. This is done by relaxing 
the particle conservation constraints on some moments; 
in our specific case, only p c and p p o are kept conserved 
while the other moments of the polymer density distri- 
bution are adapted to minimise the free energy. One 
can show that this reproduces the semi-grandcanonical 
analysis of Sear and Frenkel [24j, in full analogy to our 
discussion in Ref. 1271. 



III. PHASE BEHAVIOUR 

In this section we will describe our results for the over- 
all phase behaviour of a mixture of polydisperse poly- 
mers and monodisperse hard-sphere colloids. Our numer- 
ical work requires a choice to be made for the "parent" 
polymer diameter distribution np°^(<7 p ) which is imposed 
when the system is prepared. We concentrate primarily 
on a Schulz distribution 



n p °\a p ) oc CTpexp 
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An upper cutoff of a p = 1 on polymer sizes is imposed 
because for larger sizes the chain structure of the poly- 
mers becomes important and polymers can no longer be 
treated as effective spheres; see e.g. Refs. The 
parameters a and z control the mean and width of the 
Schulz distribution. Without a cutoff, the mean size is 
simply £s — a while the polydispersity is related to z 
via S 2 = l/(z + 1). In the presence of the cutoff these 
relations remain valid for large z, where the system is 
almost monodisperse; otherwise z and a have to be cal- 
culated numerically to give the desired values of £s and 
8. We will focus in our numerical work on three val- 
ues of the polydispersity, namely 8 = 1/^/20 + 1 ss 0.22, 
8 = 1/V5+T « 0.41 and 6 = 1/VT+T ~ 0.58. 

To check the generality of our results, calculations have 
also been performed for a triangular size distribution: 



,(()) 



_ 1 / °p — ( CT — w ) f° r a ~ w — °p — G 



(a + w) — u p for a < <j p < a 



w 



Here the mean size is £s = & independently of the width 
parameter w, which is related to the polydispersity by 
w = \f%a8. We find that qualitative features do not 
depend strongly on the shape of the polymer size distri- 
bution and so mostly report only the results for a Schulz 
distribution. 



A. Phase diagrams 

As discussed in the introduction, the phase diagrams of 
monodisperse colloid-poymer mixtures may show either 
only fluid-solid coexistence or, for sufficiently large poly- 
mers, an additional gas-liquid region with an associated 
gas-liquid-solid triangle. We begin by investigating the 
effect of polymer polydispersity on the change in phase 
diagram topology as polymer size is increased, by cal- 
culating phase diagrams for a range of values of mean 
polymer size £s and polymer polydispersity 8. Fig. [5] 
summarises the results, for Schulz size distributions. The 
empty circles indicate phase diagrams containing gas- 
liquid coexistence regions, while the full circles identify 
phase diagrams where only a fluid to fluid-solid transi- 
tion was found. From these results we can estimate the 
threshold value £ c for the mean polymer diameter as a 
function of polydispersity 8, as indicated by the thick 
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FIG. 2: Phase diagram topologies as a function of mean poly- 
mer size £s and polydispersity S. Full circles indicate param- 
eter values where only fluid-solid coexistence is found, while 
empty circles show that the phase diagram has a gas-liquid co- 
existence region. The thick solid line shows a rough estimate 
for the boundary between the two regions, which defines the 
crossover value of the mean polymer size, £c(<5). Two sample 
sets of curves corresponding to constant £s, £iv, £m an d £s 
(chosen to agree at S = 0.41) are also drawn. Comparison 
with the thick solid line shows that £m is most suitable for 
predicting the phase diagram topology independently of 5. 



solid line in Fig. [21 The data show that as polydispersity 
increases, £ c diminishes significantly. This is due to the 
fact that, for given £5, polydispersity tends to favour gas- 
liquid over fluid-solid coexistence. At £s = 0.2, for ex- 
ample, the phase diagram topology changes at 5 ~ 0.35, 
from fluid-solid coexistence only at low S to additional 
gas-liquid and gas-liquid-solid coexistence at high 5. 

In Fig. we show the full phase diagram obtained 
for Schulz distributions of polydispersity <5 = 0.41 and 
5 = 0.58 and at constant mean polymer size of £s = 0.2 
(top) and £s = 0.4 (bottom). The axis variables are the 
colloid volume fraction </> c - identical in our units to p c 
- and the (effective) polymer volume fraction tfi p — p p a, 
following the suggestion by Warren [2^ that the repre- 
sentation in terms of <p p minimises polydispersity effects 
on the phase behaviour. At = 0.2 we see explic- 
itly how gas-liquid and three-phase regions appear as S 
is increased, while the fluid-solid boundary recedes. At 
£ = 0.4, Fig. (bottom), a gas-liquid region is present for 
all polydispersities. The main effect of increasing poly- 
dispersity is again to favour gas-liquid coexistence. We 
note that, in contrast to what is found in a monodisperse 
system, the phase boundaries of the gas-liquid-solid re- 
gion are no longer straight. This is due to the fact that 
in a polydisperse system the constraint of particle con- 
servation must be satisfied for infinitely many different 
polymer sizes; indeed the polymer size distributions of 
the three coexisting phases changes as one moves around 
the three-phase region. 




FIG. 3: Phase diagrams for Schulz distributions with poly- 
dispersity S = 0.41 and 8 — 0.58 and fixed mean polymer 
size Cs = 0.2 (top) and £s = 0.4 (bottom), plotted in terms 
of effective polymer volume fraction <j> p versus colloid volume 
fraction cj> c . Dotted line: best guess for the phase boundary 
where our numerical data become unreliable. The inset shows 
an enlarged view of the three-phase region for S — 0.41 and 
Cs = 0.2. 



We next ask which measure of effective polymer size is 
most appropriate for predicting the phase diagram topol- 
ogy. Ideally, as explained in the introduction, the value of 
this quantity should be enough to predict the occurrence 
or not of gas-liquid coexistence, whatever the polydis- 
persity. In terms of Fig. |2J this means that the curve 
£ c (<5) should follow a contour of constant effective poly- 
mer size. We have therefore plotted such contours for 
each of the four possible measures of effective size dis- 
cussed in the introduction (£5: mean diameter, £j\r: num- 
ber average, £m'- mass average, £b'- virial average). Two 
sets of contours are shown, chosen such that all four quan- 
tities have identical values (0.2 and 0.4, respectively) at 
polydispersity 6 = 1/V5 + 1 ~ 0.41. This value of 5 
was also used as a reference point in |24(. We see that 
among the four candidate definitions of an effective poly- 
mer size, the mass average gives a contour which fol- 
lows most closely - although certainly not perfectly - 
the curve £ C (<S) which marks the crossover between the 
two phase diagram topologies. It is therefore our best 
candidate for a useful measure of effective polymer size. 



6 




FIG. 4: Top: Phase diagrams for Schulz distributions with 
polydispersity S = 0.22, 0.41 and 0.58 and constant mass 
average £m = 0.558; this should be contrasted with Fig. [3] 
where the mean size £s (= 0.4) was held fixed. Bottom: Phase 
diagrams for a triangular size distribution with 5 — 0.2 and 0.4 
at constant mass average £m — 0.528. In the inset the results 
for fixed mean diameter of £s = 0.4 are shown for comparison. 
Dotted lines indicate best guesses for phase boundaries we 
could not determine accurately from our numerics. 



We can now go further and investigate whether poly- 
mer size distributions with the same mass-average £m 
have not only the same phase diagram topology but in 
fact quantitatively similar phase diagrams. As pointed 
out in the introduction, the answer to this question also 
depends on how we represent the concentration of poly- 
mers; we continue to use 4> p . With this choice, the col- 
lapse of phase diagrams for polymer size distributions 
with the same £m is in fact rather good, as shown in Fig.|3] 
(top). The phase diagrams displayed correspond to the 
parameter values at the points (a)-(c) in Fig. [3 In Fig.0] 
(bottom) we show analogous results for triangular size 
distributions. Again the comparison at constant mass- 
average £m gives a good collapse of the phase diagrams 
for different polydispersities, while at fixed £s polydis- 
persity causes significant changes; see the inset of Fig. 0] 
(bottom). In summary, we can say that systems with 
the same mass-average £m give phase diagrams which are 
not only qualitatively but in fact quantitatively largely 
independent of polydispersity. Our results suggest that 



this conclusion holds independently of the polymer size 
distribution - though multimodal distributions might be 
expected to behave differently - and up to fairly substan- 
tial polydispersities of at least 40% and possibly higher. 

We note finally that the results in Fig. are for the 
case of relatively large polymers, where phase separa- 
tion is initially into gas and liquid except at very low 
polymer concentrations. We have checked the near S- 
independence of phase diagrams with the same £m also 
for smaller polymers, where fluid-solid phase separation 
occurs first (data not shown), and reached the same con- 
clusion. Very close to the crossover in the phase diagram 
topology, deviations do of course occur. This is clear 
from Fig. |2 which shows that the contour of constant £m 
does not follow the crossover curve perfectly. In a narrow 
range of fixed values of £m the phase diagram will there- 
fore change with 8, both in terms of its topology and the 
quantitative location of the phase boundaries. 



B. Comparison with semi-grandcanonical approach 

The qualitative trend of the results presented so far is 
that, for fixed mean polymer size £5, polymer polydisper- 
sity favours gas-liquid coexistence over fluid-solid phase 
splits. On the other hand, the results obtained by Sear 
and Frenkel in the semi-grandcanonical approach |2J| , as 
sketched in Fig. ^ suggested that fluid-solid phase sepa- 
ration can be favoured by polydispersity in some regions 
of the phase diagram. This difference clearly needs to be 
understood. 

Three possible causes suggest themselves. First, the 
results of Sear and Frenkel were obtained by comparing 
reservoir polymer size distributions with constant virial 
average size £b, whereas we initially focused on con- 
stant £g. Second, their phase diagrams are represented 
in terms of polymer density p p rather than effective poly- 
mer volume fraction <p p . Third, their polymer size dis- 
tribution in the system varies because it is imposed only 
indirectly via the reservoir size distribution, whereas it is 
fixed in our analysis. 

It is intuitively clear from the discussion above that 
keeping £g rather than £g constant should not influence 
the results greatly. Indeed, if we look at Fig.|21we see that 
contours of constant £b are close to those for constant £g, 
which arc vertical lines. To rule out more quantitatively 
the first two possible reasons discussed above, we have 
calculated the cloud curves for the onset of gas- liquid co- 
existence for three Schulz size distributions with identical 
£b ; and plotted them in Fig. [S] (top) with p p on the ver- 
tical axis. As before, gas-liquid coexistence is favoured 
by increased polymer polydispersity. This leaves only the 
third reason, i.e. the difference between the canonical and 
semi-grandcanonical scenarios. Fig. [S] (bottom) compares 
our results for the semi-grandcanonical scenario with the 
canonical ones in the top graph. The same three poly- 
mer size distributions are used as before, but they now 
specify the properties of the reservoir: the normalised 
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FIG. 5: Cloud curves showing the onset of phase coexistence 
coming from low colloid density. Top: Canonical scenario; 
results are shown for three parent size distributions fl p (<r p ) 
with different polydispersities, all with the same virial average 
polymer size £b as the reference distribution with = 0.4 
and 5 = 0.4L Bottom: Results for the corresponding semi- 
grandcanonical scenario, where the same three distributions 
are used as the reservoir size distribution, n P (cr p ). 



reservoir size distribution, np(<r p ), is the same as the 

parent size distribution rt p °^(cr p ) used in the canonical 
case. As expected, our semi-grandcanonical results are 
in good agreement 0] with those of Ref. [24[. In par- 
ticular, we see in Fig. [5] (bottom) that the change to the 
semi-grandcanonical description has reversed the order of 
the cloud curves at low polymer density as compared to 
the canonical scenario shown in Fig. [S] (top). 

We can thus conclude that the opposite trends with 
polydispersity seen in the present study and in Ref. 
arise from the choice of a semi-grandcanonical scenario 
in the latter. To understand explicitly how the polymer 
size distribution is affected by this, we can combine l|4*|l 
and © to write 

PpK) = PpK) ex P ( - pSm^pJ ( 9 ) 

where p p (cr p ) is the polymer density distribution in the 
reservoir as before. The excess chemical potentials /ig j 
that appear here vary as we change the colloid density 



of our system. Correspondingly the polymer density 
distribution has its shape modified as we move around 
the phase diagram. In particular, at the onset of phase 
coexistence the normalised size distribution n p (cr p ) will 
differ from the one in the reservoir, rt p (cr p ). In the 
more realistic canonical description, on the other hand, 

n p( a p) = n P°\ (T p) remains fixed up to the point where 
phase coexistence begins. 

To illustrate this difference, we plot in Fig. [fj] (top) 
the normalised size distributions n p (a p ) in the semi- 
grandcanonical system (solid lines). These are compared 
to the reservoir size distributions nl(ap) (dashed lines), 
which are the same as used in Fig. ]5\ (bottom); the col- 
loid and polymer densities are (<f> c , pp) = (0.3,1.0). We 
notice that the system and reservoir size distributions 
differ significantly: the Gibbs-Boltzmann factor in © 
shifts the mean polymer size to lower values, and this 
effect becomes stronger as the polydispersity increases. 
Fig. (bottom) traces the mean and polydispersity of 
the system size distribution as a function of colloid vol- 
ume fraction (f> c , at fixed p p = 1.0. The reduction in 
the mean size compared to the reservoir is seen to in- 
crease with colloid density. The polydispersity, on the 
other hand, is only weakly affected, as the inset of Fig- El 
(bottom) shows. 

The lines in Fig. [S] (bottom) are drawn up to the cloud 
point where phase separation begins. For the widest 
reservoir distribution (<5 = 0.58), the mean polymer size 
in the system at this point is roughly half that in the 
reservoir. Because smaller polymers produce a depletion 
interaction that is both weaker and more short-ranged, 
this reduction in polymer size delays the onset of phase 
separation. Overall, we conclude that the inversion in 
the order of the cloud curves in the bottom part of Fig. [3] 
(bottom) results from the shift of the polymer size distri- 
bution to lower values as the polydispersity of the reser- 
voir is increased. 

The reduction of the polymer mean size in the semi- 
grandcanonical scenario can be understood quantita- 
tively for narrow reservoir distributions. Looking at 10, 
if pp(crp) is narrowly peaked around a p = then we 
can expand the exponential as const X [1 — A(cr p — ££)] to 
leading order. From this one easily sees that the mean 
size £s in the system is 

£ S = e s a<k - ^) 2 ) r =e s - \s 2 (e s ) 2 (10) 

to lowest order in 5. The coefficient A must be positive 
because the excess chemical potential in a hard-sphere 
system, Pi^(c p ), is an increasing function of <j p . As ex- 
pected we have, therefore, a reduction in the mean poly- 
mer size compared to the reservoir, by an amount which 
grows with the reservoir polydispersity as 5 2 . 

Finally, we note that the three curves in Fig. [3] (bot- 
tom) all appear to cross at the same point. To see how 
this effect arises, we focus again on narrow reservoir dis- 
tributions. The location of the cloud point value of the 
polymer density p p = ,o p (0 c , £s, <5) depends on the col- 
loid volume fraction </> c , the mean polymer size £s in 
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FIG. 6: Top: Comparison of polymer size distributions in 
a semi-grandcanonical system, n p (er p ) (solid lines), with the 
corresponding reservoir size distributions np(<r p ) (dashed). 
The same three reservoir distributions as in Fig. (bot- 
tom) are used; the colloid and polymer densities are fixed 
to (0c,Pp) = (0.3, 1.0). Bottom: Mean polymer size £s (main 
graph) and polydispersity S (inset) of the system size dis- 
tribution n p (o"p), at fixed polymer density p p = 1.0 and for 
varying colloid density <f> c . Arrows on the vertical axis in the 
main plot indicate the mean sizes in the reservoir. 



the system and the polydispersity 5; for small 5, other 
properties of the size distribution are irrelevant. In a 
semi-grandcanonical setting, £s and S are determined in- 
directly by the colloid volume fraction cf> c as well as the 
mean polymer size £!j and the polydispersity of the reser- 
voir. In a perturbation expansion to order S 2 we can ap- 
proximate S 2 by its fixed value in the reservoir, and write 
£s = S,s ~ ^ 2 {Cs) 2 irom GH- The explicit dependence 
on S of the cloudpoint is also expected to be quadratic 
to leading order ^i|[5(|. Overall, we can therefore ex- 
pand the cloud point position in a weakly polydisperse 
semi-grandcanonical system as 

p P (<f> c ,£s,5) = ft.focG.O) - A<S 2 (&) 2 ^(<k,&,0) 

dt,s 

+ ^W> c ,^,o) (ii) 

How this varies with increasing polydispersity S 2 depends 
on the competition of the last two terms. The last term 
is also present in the canonical scenario and represents 



the shift in the cloud point at constant mean polymer 
size. We saw above that polydispersity favours phase co- 
existence under these conditions, so that dp p /dS 2 must 
be negative; see also Fig. [S] (top). The second term is 
only present for the semi-grandcanonical case; because 
increasing £s again favours phase separation, we have 
that also dp-p/d^s is negative. Thus the two <5 2 -terms 
in l|ll|) have opposite signs. The coefficient A depends on 
4> c ; it tends to zero for <p c — » 0, where the excess chemi- 
cal potentials in a hard-sphere system vanish, and grows 
from there with increasing <p c . At small <p c we there- 
fore expect the third term in 111|) to dominate, so that 
polydispersity favours phase separation as in the canon- 
ical case. As <fi c increases, the second term eventually 
balances the third. At this point the cloud curve is un- 
affected by 6 2 to leading order, and this produces the 
crossing phenomenon in Fig. [3] (bottom). At even higher 
(j) c , the second term dominates and so phase separation 
is delayed with increasing 5, exactly as seen in the lower 
part of Fig. (bottom). 

Overall, while in the realistic canonical description 
polydispersity always favours gas-liquid phase separa- 
tion, in the semi-grandcanonical scenario this trend can 
be counteracted by a decrease in the mean polymer size 
in the system. This effect can become dominant at high 
colloid volume fractions, producing the reversal in the 
order of the curves in the lower parts of Fig. [3] (top) and 
(bottom). 



IV. CONCLUSION AND OUTLOOK 

We have investigated the effect of polymer polydisper- 
sity on a mixture of ideal polymers and hard-sphere col- 
loids. Our focus was on the realistic (canonical) case 
where the distribution of polymer sizes is fixed when the 
system is prepared, but we have also compared with semi- 
grandcanonical approaches where the size distribution 
adapts by equilibrating to chemical potentials imposed 
by a large polymer reservoir p^j . For the canonical set- 
ting, gas-liquid phase separation is favoured by increased 
polydispersity while fluid-solid phase separation is re- 
tarded. This is in broad agreement with the intuition 
that polydispersity disfavours ordered phases 0, |27| . 
though less obvious here since the size polydispersity is 
in the polymers, not in the colloids that order transla- 
tionally in the solid. 

In a semi-grandcanonical scenario the above trends 
with polydispersity are reversed in some parts of the 
phase diagram 24]. We saw that this arises because 
the polymer size distribution in the system can become 
rather different from the reservoir size distribution, es- 
pecially for dense colloidal phases. Specifically, the size 
distribution is shifted to smaller polymer sizes, the more 
so the larger the degree of polydispersity. This effect is 
sufficiently strong to reverse the order of the cloud curves 
compared to the canonical setting, favouring fluid-solid 
over gas-liquid formation. 
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For the canonical case of conserved polymer size distri- 
bution, we established that the phase diagram topology 
was nearly independent of polydispersity provided one 
compares systems with the same mass average polymer 
size £m- We were lead to this observation by estimating 
the location of the curve £, C (S) in the (£s, 8) plane along 
which the phase diagram topology changes, and com- 
paring this with contours of constant £m- The crossover 
value £ c decreases with increasing polydispersity <5, show- 
ing quantitatively that gas-liquid coexistence is favoured 
as polydispersity is increased at constant mean polymer 
size. Finally, we showed that even quantitatively the 
phase diagrams of systems with the same mass average 
£m are nearly independent of polydispersity, provided 
that polymer concentration is represented in terms of the 
effective polymer volume fraction 4> p . 

To assess the reliability of our results, we discuss briefly 
the approximations made. We have employed the stan- 
dard AO model instead of explicitly modelling ideal poly- 
mer chains. Simulations suggest that this is a rea- 
sonable approximation for polymers which are not too 
large, i.e. <j p < 1. A slight improvement can be obtained 
by using a mapping from polymer radius of gyration to 
effective sphere radius er p /2 which accounts for the effect 
of the curvature of the colloid surface on the polymer 
conformations [n|,[£l|. In the limit of a small polymer, 
this results in an increase of a p by a factor rs 1.13 

over the naive value, but in the range 0.2 < a p < 1 of 
interest the corrections to a p are small (< 5%). With 
such corrections included, a Schulz distribution across 
the polymers' radii of gyration - the case we concen- 
trated on - would result in a distribution of <7 P that is no 
longer precisely of Schulz form; quantitatively, however, 
this will again be a small effect. We have also used the 
free- volume approximation to estimate the free ener gy o f 
the polydisperse AO model. As discussed elsewhere [21j . 
this is essentially equivalent to the density functional the- 
ory approach of Ref . [33| ; see also the discussion for the 
monodisperse case in Ref. . Again, simulations sug- 
gest that this approach is reliable 0, |2(| in the range 
0.3 < dp < 1, although somewhat less so after polymer 
reservoir densities p p are converted back to those in the 
system (compare e.g. Fig. 3 in Ref. [23 with Fig. 2(b) in 
Ref. E3). 

We comment briefly on comparisons with experimen- 
tal work. Much work on mixtures of colloids and 0-point 
polymers has been done in Edinburgh; see e.g. the re- 
view J34(. Experimentally, the typical polymer size is 
normally determined from the weight-average molecular 
weight, which corresponds to our Interestingly, our 
results then suggest that polymer polydispersity should 
have only a small effect on the phase behaviour, so that 
e.g. the comparison between experimental data and AO 
model predictions in jE^ would remain essentially un- 
changed |5^| . A systematic experimental study of such 
polydispersity effects would seem worthwhile. For ex- 
ample, one could prepare two polydisperse systems, one 
whose number-averaged molecular weight (~ and one 



whose mass-averaged molecular weight (~ £f^) coincides 
with that of a near-monodisperse reference system; our 
theory would then predict that the second polydisperse 
system produces phase behaviour much more similar to 
the monodisperse reference than the first. Also fruit- 
ful could be the investigation of polymer size fractiona- 
tion between coexisting phases. We have not displayed 
our theoretical predictions for this; one finds as in the 
semi-grandcanonical approach [24| that the larger poly- 
mers are found in the phases with lower colloid density, 
but the detailed distribution shapes are different because 
they always need to combine to give the fixed parent dis- 
tribution. 

An interesting extension of the present work would be 
to consider the combined effect of polymer and colloid 
size polydispersity, generalising the present study and 
our previous one on polydisperse colloids and monodis- 
perse polymers plj . This is in principle possible using 
the current framework, but the excess free energy would 
then depend on eight moments p C i and p P i {i = . . .3), 
making the problem very challenging numerically. Phys- 
ically, one would expect from our current results that 
the crossover value of the polymer size, when measured 
in terms of the mass average £m> would be largely un- 
affected by polymer polydispersity, but decrease signifi- 
cantly with increasing colloid polydispersity plj . A sim- 
pler alternative to the treatment of a fully polydisperse 
system would be to consider the colloid size distribu- 
tion as "quenched" , i.e. equal to the parent in all phases, 
thus excluding colloid size fractionation. One would then 
need only one conserved colloid moment in the excess 
free energy, and the complexity of the numerical analy- 
sis would be the same as in the present study. Physi- 
cally, this quenched approximation could be appropriate 
for describing the initial stages of phase separation in a 
polydisperse colloid-polymer mixture, where polymer size 
fractionation - which should be the faster process - has 
already taken place while colloid size fractionation is still 
negligible. 

A key open question for further work is how our re- 
sults generalize to the case of good or poor (rather than 
9-) solvents, where polymer-polymer interactions can no 
longer be neglected. For monodisperse polymers this 
issue has been tackled by a variety of methods includ- 
ing virial expansions near the 0-point |54j , integral equa- 
tions [53, llg, |53, ISg , Flory-Huggins theor y 1591. coarse- 
graining of polymers into soft colloids [20I l5ll l60||. mod- 
ifications of free- volume theory [§l|, and an AO model 
with added soft polymer-polymer repulsions |62l ] to model 
good solvents or with a repulsive co-solvent to mimic 
poor ones 63]. It remains an open challenge to in- 
corporate polymer polydispersity into these approaches. 
Some progress in this direction has already been made in 
Ref. [2{| , which was published after the present work was 
completed. In this study the polymers were modelled as 
chains of hard-sphere particles, as would be appropriate 
for good solvent conditions, and calculations for the onset 
of gas-liquid coexistence were performed using a truncat- 
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able model free e nerg y derived from thermodynamic per- 
turbation theory [64( • Also in this scenario with interact- 
ing polymers, it is found that polydispersity significantly 
enlarges the size of the coexistence region. It appears 
therefore that this qualitative conclusion is robust to the 
presence or otherwise of polymer non-ideality. If the ap- 
proach of Ref. [2f| can be extended to include colloidal 



solids, it would be interesting to ask whether polydis- 
persity effects on the relative stability of fluid-solid and 
gas-liquid phase separation are likewise robust under the 
inclusion of polymer-polymer interactions. 

The authors acknowledge support of the EPSRC 
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